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INTRODUCTION 


It  is  generally  accepted  that  some  mechanism  for  ionizing  the 
neutral  gas  near  an  ionospheric  probe  at  high  positive  potential  is 
required  to  account  for  observed  levels  of  electron  current  collection 
(Winkler,  1980).  Early  theoretical  studies  (Beard  and  Johnson,  1961; 
Parker  and  Murphy,  1967;  Linson,  1969)  which  considered  only  the 
collection  of  ambient  electrons  tended  to  underestimate  the  return 
currents  to  an  electron  emitting  vehicle  or  conversely,  overestimate  the 
charge  up  potential.  The  predictions  of  these  studies  varied  widely 
depending  primarily  upon  the  inclusion  of  magnetic  field  effects,  with 
Beard  and  Johnson  ignoring  the  geomagnetic  field;  Parker  and  Murphy 
including  it,  but  mostly  ignoring  space  charge  effects,  and  Linson 
proposing  "turbulence"  to  transport  electrons  across  the  magnetic  field. 
The  ionization  of  neutrals  has  been  considered  by  Zhulin  et  al .  (1976) 
and  Galeev  et  al .  (1976)  but  the  only  attempts  at  simultaneously 
including  ionization  and  space  charge  (but  not  magnetic  fields)  are  the 
related  studies  of  Leadon  et  a  I .  (1980)  and  Lai  et  a  I .  (1983). 

The  present  study  is  based  on  the  configuration  shown  in  Figure  1, 
where  the  electrons  attracted  by  the  positive  probe  fall  through  the 
sheath  region  and  impact  ionized  neural  atoms  with  a  cross-section  a. 

The  newly  liberated  secondary  ions  move  out  through  the  sheath  and  the 
secondary  electrons  are  collected  by  the  probe.  The  secondary  charge 
densities  should  be  self-consistent  with  the  solution  of  Poisson's 
equation.  Lai  et  al .  predict  a  non-monotonic  current-voltage 
characteristic  that  seems  to  be  capable  of  explaining  many  observations, 
but  the  inclusion  of  magnetic  field  effects  may  ultimately  be  required. 


The  present  study  was  prompted  by  a  numerical  catastrophe  that 
occurs  in  the  Lai  model  whenever  the  sheath  potential  profile  becomes 
non-monotonic.  This  problem  stems  from  the  use  of  analytic  integral 
expressions  for  the  secondary  space  charge.  In  our  model,  the  secondary 
ion  space  charge  is  treated  by  a  fluid  model  that  overcomes  this 
problem.  With  our  model,  we  have  been  able  to  model  sheath  behavior  as 
the  secondary  ion  density  in  the  (otherwise  electron)  sheath  reaches  and 
exceeds  a  critical  level'.  Specifically,  the  critical  effect  is  that 
with  sufficient  internal  ion  space  charge,  the  sheath  edge  moves 
outwards  to  include  more  negative  electron  space  charge,  but  in  doing 
so,  the  interior  ion  production  is  proportionally  increased  leading  to 
continued  sheath  expansion,  or  "explosion”.  This  effect  is  similar  to 
the  propagation  of  double  layers  (Carlqvist,  1972). 


THE  SPACE  CHARGE  LIMITED  SHEATH 


This  study  of  sheath  ionization  is  focused  on  the  regime  of  space 
charge  limited  current  collection.  This  is  equivalent  to  the  short 
Debye  length  limit  <<  where  Xq  is  the  Debye  length  and  is  the 
probe  radius.  We  also  presume  that  the  probe  voltage  Vp  »  kT/e.  Under 
these  conditions  a  plasma  sheath  may  be  described  as  composed  of  three 
regions:  the  sheath,  the  sheath  edge,  and  the  pre-sheath.  The  sheath 
is  the  region  adjacent  to  the  probe  where  only  the  attracted  specie  is 
present  from  the  ambient  plasma,  and  the  potential  decreases  with 
radius,  r,  faster  than  r~'^.  This,  combined  with  the  assumption  of  high 
probe  potential,  produces  a  strong  electric  field  that  overcomes  any 
orbital  motion,  and  results  in  near  radial  trajectories  inside  of  some 
absorption  radius.  The  sheath  edge  is  the  nominal  absorption  radius  for 
the  bulk  of  the  attracted  distribution,  and  is  generally  a  few  Xq  in 
thickness.  The  presheath  is  the  quasi  neutral  region  that  connects  the 
sheath  to  the  undisturbed  plasma  (Parrot  et  a  I . .  1982). 

It  is  within  this  picture  of  the  sheath  that  we  propose  to  examine 
sheath  structure  using  space  charge  limited  diode  models  (Chen,  1965). 
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SATURATED  DOUBLE  DIODES 


The  objectives  of  this  study  were  two-fold:  one,  determine  the 
relationship  between  the  probe,  plasma,  and  neutral  background  that  give 
rise  to  an  explosive  condition;  and  two,  model  the  evolution  of  the 
sheath  when  parameters  have  exceeded  the  threshold  for  sheath  explosion. 
The  emphasis  is  on  the  first  objective.  Thus  we  hypothesize  that  sheath 
ionization  is  simulated  by  planar  and  spherical  double  diodes,  then 
verify  and  quantify  this  assumption  with  our  numerical  model.  By  double 
diode  we  mean  a  two  electrode  device  where  charge  carriers  of  opposite 
sign  are  emitted  from  the  opposing  electrodes,  and  transported  across 
the  gap  between  the  electrodes. 

The  one  double  diode  model  to  which  solutions  can  be  found  in  the 
literature  is  due  to  Langmuir  (1929),  and  generalized  to  double  layers 
by  Block  (1972).  This  model  is  i 1 1 ustrated ' i n  Figure  2  as  a  double 
layer  with  anode  and  cathode  locations  indicated  to  make  the  connection 
with  Langmuir’s  diode.  We  assume  a  strong  layer  where  the  initial  or 
thermal  velocities  are  negligible  compared  to  potential  drop.  If  the 
anode  does  not  emit,  then  we  have  the  usual  planar  Chi Id-Langmui r  diode 
(Child,  1911;  Langmuir,  1913)  where  the  separation  D^|_,  potential  drop 
V,  and  current  density  J  are  related  by  the  expression 


J 

e 


or,  if  we  use  for  Jg,  the  plasma  thermal  current, 


Je  =  Ne('<T/2)rmg) 


we  may  wr  i  te 


Dcl  =1.26  XD(eV/kT)3/4 


(2) 


(3) 
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where  Xq  s  is  the  plasma  Debye  length.  If  we  allow  a  space 

charge  limited  ion  current  to  flow  from  the  anode,  we  have  a  symmetric 
potential  profile,  and  for  the  same  electron  current,  the  separation 
distance  increases  to 

D<i  =  1.36  D^l, 

or  if  the  diode  gap  were  fixed,  the  electron  and  ion  currents  would  be 
1.86  times  the  limit  (2).  When  both  specie  currents  are  space  charge 
limited  we  see  from  (2)  (regardless  of  the  value  of  D)  that 

Je  =  '''"i/'"e  Ji  (*) 

This  relation,  known  as  the  "Langmuir  condition",  is  the  basic  stability 
requirement  for  a  strong  double  layer.  If  either  flux  should  deviate, 
one  sign  of  space  charge  in  the  layer  will  dominate,  and  the  layer  will 
move  upstream  of  the  weaker  flux  so  as  to  satisfy  this  condition  in  the 
moving  frame  of  reference. 

This  is  all  well  and  fine,  but  in  our  picture  of  sheath 
ionization,  the  ions  are  produced  throughout  the  sheath  region  and 
rather  than  the  probe  surface.  To  our  knowledge  the  literature  does  not 
contain  a  complete  solution  for  such  a  diode,  although  Langmuir  (1929) 
did  analyze  the  problem  for  single  ion  release. 

The  situation  for  a  diode  formed  of  concentric  spherical 
electrodes  is  even  more  barren.  We  do  have  available  the  work  of 
Langmuir  and  Blodgett  (1924)  who  studied  the  single  emitting  electrode 
spherical  diode.  Their  results  are  primarily  tabular  and  are  not 
reproduced  here,  but  we  note  that  Parker  (1980)  provides  a  useful  fit  to 
those  results.  An  interesting  feature  of  their  re.'<ults  is  that  for  a 
given  potential  drop  and  sphere  radii,  the  limiting  current  (electrons) 
is  higher  for  a  central  cathode  than  for  a  central  anode.  This  is  not 
surprising  since  mutually  repulsive  particles  are  happier  diverging  than 
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converging.  This  does  lead  us  to  suspect  that  the  double  spherical 
diode  will  not  have  the  symmetry  that  allowed  us  to  jump  from  Eq.  (2)  to 
Eq.  (4),  and  that  a  flux  of  ions  flowing  from  an  inner  anode  will  exceed 
that  suggested  by  (4)  by  an  amount  that  we  call  the  anode  factor  a,  thus 

F.  =  aVm  /m.  F  , 

I  e'  I  e  ' 

where  a  wi 1 1  be  determined  numerically. 


EXPLOSIVE  SHEATH  IONIZATION 

The  concept  of  explosive  sheath  ionization  can  be  illustrated 
quite  simply.  Consider  a  spherical  probe  of  radius  r^,  with  a  sheath  of 
volume  A,  area  A^,  and  an  average  ion  production  rate 


ff(6(r))  F^(r)  •  Airr^  dr 


where  r^  is  the  radius  at  which  electrons  have  attained  sufficient 
energy,  £,  for  the  ionization  of  neutrals  with  density  N^.  Assuming 
that  orbital  motion  within  the  sheath  is  negligible,  we  may  write 


F  (r)  =  F  f^1 
e'  ‘  es  Ir  J 

where  F^^  is  the  presheath  enhanced  (Parrot,  1980)  ambient  electron 
flux  to  the  sheath  edge.  By  also  assuming  that  is  constant  with  e, 

we  have 


5(rt)  =  ^rrg^  A"^  F^gCr^  -  r^,)  . 


The  ion  flux  out  through  the  sheath  will  be  approximately 


F  =  S(rJ  T-  =  aN  F  (r  -  r  ) 
I  '  t''  Ag  n  es'  t  o' 


We  claim  that  if 


F.  >  aJm/m. 


then  the  sheath  will  move  outwards,  >0,  r^  >  0,  leading  directly  to 
Pj  >  0,  resulting  in  continued  instability  of  the  sheath.  From  Eq.  (5) 
we  can  see  that  condition  6  is  equivalent  to 


where  Xjg  =  is  the  mean  free  path  of  electrons  for  ionization  of 

neutrals  and  D3  -  (r^  -  r^)  is  the  sheath  thickness.  A  careful  estimate 
of  Tg  is  difficult  to  obtain  since  the  ion  distribution  function  will 
depend  on  the  specific  potential  distribution  within  the  sheath,  and  the 
energy  dependence  of  the  ionization  cross-section.  However,  we  may 
typify  the  expansion  by  noting  tl  ^t  it  will  be  limited  by  the  speed  of 
fastest  ions,  so 


"■s 


< 


2e  AV 
m 


m/s 


for  oxygen  atoms  and  a  probe  potential  of  50  volts. 


There  should  certainly  be  mechanisms  that  will  modify  and  limit 
the  simple  analysis  present  here.  Some  of  these  might  be: 


A  complex  sheath  potential  profile  may  modify  ion  velocity 
distributions  and  cause  some  to  impact  the  probe.  Much  of  this  effect 
will  be  included  in  the  following  numerical  model. 

Anomalous  resistivity,  which  if  we  do  not  suggest  it,  someone  else 
will,  and  could  retard  electrons  in  the  sheath  and  increase  their  space 
charge . 


Trapped  orbits,  which  might  be  quickly  filled  by  scattered 
electrons.  These  will  not  be  truly  stable  orbits  unless  there  exists 


ii 

I 

i'.  • 


It 


f 


I  -r*.  --fc  '.-v 


%  S  %' 


regions  where  the  potential  falls  off  slower  than  r~  (Laf ramboi se, 

1966)  but  our  numerical  studies  indicate  that  as  a  sheath  becomes 
explosive  such  regions  may  indeed  appear,  even  though  the  problem  was 
initially  posed  in  the  strong  space  charge,  short  Debye  length  limit. 

Neutral  depletion  should  have  a  minimal  effect.  The  depletion 
rate  is  -  5  -  ?  ^e^s  “  (^0"^  *"*  10*^  sec“^)N5  for  a  =  2x10"^®  m^,  and 

a  range  of  possible  ionospheric  electron  fluxes,  and  where  is  the 
neutral  density  in  the  sheath.  The  repletion  rate  should  be 
approximately  ftp  -  V^^N^/r^  -  (10^  *■*  10®  sec”^)Np  for  1-10  m  sheaths, 
and  100-1000  *K  neutral  temperature  range,  so  we  see  that  Ng  =  Np. 


I 

r.  ■ 


Probe  impedance  will  have  a  significant  effect  and  is  probably 
the  major  distinction  between  this  study  and  that  of  Lai  et  a  I .  They 
modeled  their  probe  as  part  of  constant  current  circuit  where  the  probe 
voltage  would  fall  as  a  result  of  the  net  increase  in  electron  current 
that  a  sheath  explosion  would  produce.  Our  probe  is  a  constant  voltage, 
zero  impedance  device  which  can  support  arbitrary  large  currents,  which 
is  a  bit  removed  from  the  case  of  an  electron  emitting  rocket,  but  does 
relate  to  a  small  probe  on  a  large  vehicle  biased  by  a  hefty  power 
supply . 


From  an  electric  circuit  viewpoint,  this  sheath  explosion  would 
appear  as  a  collapse  of  the  sheath  impedance. 


IT 


r 

i 
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THE  NUMERICAL  MODEL 


A  one-dimensional  fluid  Poisson  code,  FLOMO,  was  written  to  model 
the  sheath  ionization  problem.  FLOMO  can  model  either  a  single  planar 
coordinate  or  a  spherical  radial  coordinate.  A  number  of  gridding 
schemes  were  considered,  but  uniform  x  (planar)  and  r  (spherical) 
spacings  proved  to  be  the  most  practical.  The  Poisson-sol ver  is  a 
simple  tri -diagonal  method  (Richtmyer  and  Morton,  1967)  that 
incorporates  the  following  charge  densities: 


The  repel  led  ambient  ion  density  was  modeled  everywhere  as 
N3I  =  Ng  exp(-e  V/kT) .  For  the  calculation  of  the  attracted  ambient 
electron  density,  a  sharp  edge  sheath  approximation  was  used  to  separate 
different  approximations  that  are  matched  at  the  sheath  edge.  For 
spherical  geometry  the  sheath  edge  conditions  were  taken  from  Parrot  et 
a  I .  (1982)  who  found  that  for  Vp  >>  kT/e,  the  sheath  edge  had  a 
potential  Vg  =  0.49  kT/e,  number  density  N^g  *  0.611  N^,  and  flux 
density  F^g  =  1.45  F^^.  Outside  of  the  sheath  the  attracted  electron 
number  density  was  taken  from  a  geometric  fit  to  the  results  of  Parrot 
et  a  I .  We  start  with 


as  the  geometric  density  that  would  result  from  shadowing  by  a  sheath 
subtending  a  solid  angle 


n(r)  =  If 


-  ■  1  - 


(Rg/'‘)‘ 


At  the  sheath  edge  Ngg  =  0.5  N^.  A  good  fit  to  the  results  of  Parrot  et 
a  I .  is  given  by 


l' 


where  F  s  qE/m  is  the  single  particle  acceleration  and  E  =  VV.  The 


closure  of  Eqs.  (9)  and  (10)  is  provided  by  Poisson’s  equation.  These 
equations  are  solved  by  finite  differences  using  full  upwind  and 
backwards  time  centering  (explicit)  on  the  convective  terms  (Richtmeyer 
and  Morton,  1967).  This  approach  was  chosen  because  of  its  robustness 
and  its  ability  to  handle  non-monotonic  potential  distributions,  which 
is  the  particular  advantage  over  an  integral  kinetic  method.  More 
rigorous  methods  such  as  a  turning  point  formalism  (Laf ramboise,  1966) 
or  a  fully  differenced  Vlasov  approach  (Rotenberg,  1983)  was  beyond  the 
scope  of  this  work.  We  have  performed  an  analysis  to  determine  the 
extent  of  error  introduced  by  our  fluid  model. 


Consider  a  radial  probe  of  negligible  radius  and  a  potential 
distribution  that  is  linearly  decreasing  with  r,  such  that  the  force  is 
constant  with  r.  If  we  add  to  this  a  constant  source  rate,  then  Eqs. 
(9)  and  (10)  may  be  solved  to  give 


N  s 


7 

18 


—  •  S 


f  a  0.62  S  ^ 


and 


u 


•  JFt  =  0.53  JTr 


Kinetical ly  we  have 


N(r) 


r 

f  r 

f(r,u)du  = 

S* 

r_ 

Jo  J 

0 

r 

2 


dr ' 

‘i(r,r') 


where  u(r,  r')  =  j2F(r-r')  .  This  integrates  to 


N,^  =  0.75  S  p 


12 


and  similarly  we  compute 


f(r,u)du  =0.44  JFt 


We  conclude  from  this  simple  analysis  that  the  fluid  approach  typically 
introduces  an  acceptable  error,  and  we  have  the  reassuring  result  that 


Nu  =  N|j  u. 


The  secondary  electron  space  charge  is  ignored  unless  local 
trapping  occurs.  We  justify  this  by  observing  that  in  the  absence  of 
trapping,  the  electron  space  charge  is  less  than  the  ion  space  charge  by 
a  factor  of  Vm^/mj .  Throughout  the  iterative  solution  of  our  equations 
we  inspect  the  potential  distribution  for  the  appearance  of  an  electron 
trap,  or  well.  Let  us  assume  that  a  well  of  depth  H  hns  appeared,  and 
that  in  the  well  the  electron  distribution  is  Maxwellian  for  all  trapped 
electron  energies,  zero  for  untrapped  energies,  and  normalized  to  the 
central  ion  density,  NsjC^w)*  t^'us 


N,.{r) 


N  .(r  ) 

SI  '  w' 


m 


2irkT 


3/2 


exp(e(V(r)-V(r 


))AT) 


•  4t 


-mu^/2kT  2. 
e  '  u  du 


=  N  .(rj  .  exp((V(r)  -  V(rJ)/kT 


A  secondary  electron  temperature  of  1  eV  was  assumed  in  the 
calculations.  These  various  charge  densities  are  summed  and  possibly 
adjusted  according  to  a  charge  stabilization  algorithm  (Cooke,  1983) 
that  guarantees  stability  in  the  Poisson  solution  without  iterate  mixing 
(Parker  and  Sullivan,  1969,  1970).  Briefly,  this  algorithm  derives  from 
a  plausibility  argument  which  recognizes  that  in  the  discretization  of 
Poisson's  equation,  space  charge  features  (such  as  the  sheath  edge)  that 
are  nonlinear  with  position  and  potential,  and  numerical  errors  may  both 
become  artificially  amplified  when  summed  over  discrete  elements.  The 
charge  stabilization  technique  inspects  for  excessive  total  discrete 
charge  and  reduces  it  to  a  limit  that  is  consistent  with  numerical 
stabi I i ty . 

Our  iteration  scheme  was  to  simply  timestep  the  fluid  equations 
(9)  and  (10)  for  generally  a  few  milliseconds,  iteratively  solve 
Poisson’s  equation,  and  return  to  the  fluid  equations. 


RESULTS 


The  present  results  consist  of  six  models:  A  planar  configuration 
with  either  no  local  (secondary)  ions,  an  anode  ion  source,  or  volume 
ion  production  (ionization  of  neutrals);  and  a  spherical  configuration 
with  the  same  local  ion  options.  The  planar  model  with  anode  ions  and 
no  local  ions  serves  as  a  check  on  our  numerics  since  those  results  may 
be  compared  with  the  analytic  theory  of  Langmuir  (1929).  These  models 
generally  employed  100  grid  points  to  span  30  meters  (planar)  or  12 
meters  (spherical).  All  six  models  employed  the  same  ambient  plasma 
parameters  which  are  summarized  in  Table  1.  The  steady  state  results 
are  presented  in  Table  2,  where  the  sheath  edge,  Dg,  is  found  from  the 
solution  as  the  ~0.5  kT/e  point.  The  results  for  both  the  anode  plasma 
and  ionization  models  represent  the  saturation  limit,  i.e.,  if  in  either 
case  the  ion  production  is  increased,  the  sheath  edge  begins  to  move 
outwards.  The  same  is  true  for  the  spherical  models.  The  sheath 
thickness  given  for  the  spherical  models  is  D3  s  R3  -  Rq.  The 
theoretical  prediction  for  O3  comes  from  Langmuir  and  Blodgett  (1924) 
with  the  presheath  current  enhancement  of  Parrot  et  al .  (1982). 

The  explosive  character  of  the  sheath  was  investigated  for  the 
four  cases  that  involved  saturation  limits,  but  we  present  here  a  study 
of  only  the  spherical  model  with  critical  ionization.  Figure  3  is  a 
voltage  versus  radius  profile  for  the  2  meter  radius  spherical  probe, 
with  a  neutral  density  of  5.0  x  10  .  Notice  that  the  sheath  radius  is 
>6.4  meters  which  is  essentially  the  results  for  no  ionization.  The 
charge  densities  are  presented  in  Figure  4,  where  the  attracted  electron 
density  may  be  seen  to  dominate  everywhere,  except  just  outside  the 
sheath.  This  artifact  is  primarily  a  numerical  error  relating  to  the 
precision  with  which  the  sheath  edge  may  be  located  for  the  fit  to  the 
Parrot  et  al ,  results.  This  may  also  be  observed  in  the  filtered  total 
charge  density  (labeled  as  QUSD  for  "Q  that  was  used")  shown  in  Figure 
5.  The  QUSD  values  are  normalized  to  the  ambient  plasma  density.  From 
a  comparison  of  Figures  5  and  6,  one  can  see  how  the  positive  overshoot 
error  is  limited  by  the  charge  stabilization  process.  The  recognition 
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of  this  overshoot  as  an  error  is  accomplished  by  comparison  with  the 
potential,  and  the  enforcement  of  a  charge  limit  that  is  linear  with 
potential . 

In  Figure  6,  the  neutral  density  has  been  increased  to 
Np  =  7.0  X  10*''.  The  increase  is  modest,  but  the  effect  is  dramatic 
resulting  in  an  increase  in  Rg  from  6.4  to  7.7  meters,  with  a  45  percent 
increase  in  sheath  area  and  ambient  electron  current.  Figure  7  shows 
the  charge  density  profiles  where  we  see  that  electrons  are  still 
dominant  inside  the  sheath. 

If  the  neutral  density  is  incremented  to  8.0  x  10^^  m~^,  the 
sheath  explosion  begins.  Figure  8  is  a  voltage  profile  snapshot,  at  1.3 
msec,  and  Figure  9  is  the  charge  profile.  Notice  how  the  ion  density 
begins  to  dominate  in  some  regions.  As  this  occurs,  the  electric  field 
that  removes  the  ions  is  weakened,  ion  removal  is  reduced,  and  the 
density  increases.  The  calculation  was  halted  at  2.2  msec  when  the 
sheath  hit  the  outer  boundary  at  14  meters.  This  final  set  of  profiles 
are  shown  in  Figures  10,  11  and  12. 
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CONCLUSIONS 

This  study  of  sheath  ionization  has  been  composed  of  many 
approximations  which  would  allow  one  to  question  the  accuracy  of  the 
results.  However,  our  ability  to  match  with  other  methods  suggest  that 
the  accuracy  is  adequate.  Furthermore,  while  the  determination  of  the 
critical  parameters  is  subject  to  our  numerical  accuracy,  the  basic 
instability  is  not,  and  we  are  led  to  the  conclusion  that  with 
sufficiently  low  probe  impedance,  and  sufficiently  high  neutral  density, 
an  electron  collecting  sheath  will  become  explosive.  We  may  also 
conclude  that  for  a  neutral  density  that  is  a  factor  of  2  or  so  below 
the  critical  level,  ionization  has  little  impact  on  electron  collection. 

The  determination  of  the  equivalent  anode  factor  for  the  case  of  a 
spherical  probe  with  sheath  ionization  is  subject  to  our  numerical 
accuracy,  and  this  interesting  number  deserves  closer  scrutiny. 

However,  we  may  still  obtain  a  conservative  estimate  of  the  critical 
neutral  density  from  the  planar  results  since  the  anode  factor  is  unity 
(the  limit  of  «  Rq) ,  and  the  planar  sheath  is  thickest.  Thus  we 
return  to  Eq.  (4)  and  (5),  using  S  =  o  N„  and  A/A3  -  Dj,  to  obtain 


N„  i  (.  0,)-> 


or  us i ng  Eq .  (3) , 


N„  1  {.-X,,)-* 


leVj  Jm. 


(11) 
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For  a  daytime  shuttle  environment,  we  might  choose  kT/e  s  0.1  eV, 

Ng  =  1  X  10^^  m“^,  Xq  =  0.74  cm,  and  a  =  2  x  10"^®  m"^.  For  these 
conditions  Eq.  (11)  predicts  that  a  surface  at  a  potential  of  1  kilovolt 
would  develop  an  explosive  sheath  at  ^  2  x  10^®  m“^  2  5  x  10"^  Torr. 
In  both  this  case,  and  the  previous  model  plasma,  the  critical  neutral 
density  is  quite  high  compared  to  the  natural  background,  but  within  the 
possible  contamination  range  for  the  shuttle. 
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Figure  9.  Normalized  charge  density  profiles  corresponding  to  Figure 
8.  Note  that  the  ion  density  exceeds  that  of  the  electrons 
near  mid-sheath. 
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Figure  10.  Potential  profile  at  the  conclusion  of  the  calculations  at 
2.2  msec.  The  sheath  has  expanded  out  to  the  model  boundary 
at  14  meters. 
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